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Abstract — In this paper, an idea of the cutting plane method 
is employed to improve the fractional distance of a given binary 
parity check matrix. The fractional distance is the minimum 
weight (with respect to £i -distance) of vertices of the fundamental 
polytope. The cutting polytope is defined based on redundant 
rows of the parity check matrix and it plays a key role to eliminate 
unnecessary fractional vertices in the fundamental polytope. We 
propose a greedy algorithm and its efficient implementation for 
improving the fractional distance based on the cutting plane 
method. 

I. Introduction 

Linear programming (LP) decoding proposed by Feldman 
[3] is one of the promising decoding algorithms for low 
density-parity check (LDPC) codes. The invention of LP 
decoding opened a new research field of decoding algorithms 
for binary linear codes. Recently, a number of studies on LP 
decoding have been made such as [13][7][8]. 

LP decoding has several virtues that belief propagation (BP) 
decoding does not possess. One of the advantages of LP 
decoding is that the behavior of the LP decoder can be clearly 
understood from a viewpoint of optimization. The decoding 
process of an LP decoder is just a minimization process of 
a linear function subject to linear inequalities corresponding 
to the parity check conditions. The feasible set of this linear 
programming problem is called the fundamental polytope [3]. 
The fundamental polytope is a relaxed polytope that includes 
the convex hull of all the codewords of a binary linear code. 
Decoding performance of LP decoding is thus closely related 
to geometrical properties of the fundamental polytope. 

Another advantage of LP decoding is that its decoding per- 
formance can be improved by including additional constraints 
in the original LP problem. The additional constraints tighten 
the relaxation and they lead to improved decoding perfor- 
mance. Of course, additional constraints increase the decoding 
complexity of LP decoding but we can obtain flexibility to 
choose a trade-off between performance and complexity. The 
goal of the paper is to design an efficient method to find 
additional constraints that improve this trade-off for a given 
parity check matrix. 

In this paper, we will propose a new greedy type algorithm 
to enhance the fractional distance based on the cutting poly- 
tope based on the redundant rows that eliminates unnecessary 
fractional vertices in the fundamental polytope. The additional 



constraints generated by the proposed method is based on 
redundant rows of a parity check matrix. In other words, 
the proposed method can be considered as a method to 
generate redundant parity check matrices from a given original 
parity check matrix. Therefore, the present work has close 
relationship to the works on elimination of stopping sets (SS) 
by using redundant parity check matrices and on stopping 
redundancy such as [12] [1] [6]. 

The cutting plane method is a well established technique for 
solving an integer linear programming (ILP) problem based 
on LP [10]. The basic idea of the cutting plane method is 
simple. In the first phase, an ILP problem is relaxed to an LP 
problem and then it is solved by an LP solver. If we get a 
fractional solution (i.e., a vector with elements of fractional 
number), a cutting plane (i.e., an additional linear constraint) 
matched to the fractional solution is added to the LP problem 
in the second phase. The cutting plane is actually a half space 
which excludes the fractional solution but it includes all the 
ILP solutions. In the third phase, the extended LP problem 
is solved and the above process continues until we obtain an 
integral solution. 

The fractional distance of a binary linear code is the £i- 
weight of the minimum weight vertex of the fundamental 
polytope. The fractional distance is known to be an appropriate 
geometrical property which indicates the decoding perfor- 
mance of LP decoding for the binary symmetric channel 
(BSC). It is proved that the LP decoder can successfully 
correct bit flip errors if the number of errors is less than 
half of the fractional distance [3]. In contrast to the minimum 
distance of a binary linear code, we can evaluate the fractional 
distance efficiently with an LP solver. Efficient evaluation of 
the fractional distance is especially important in the proposed 
method. 

The idea for improving LP decoding performance using 
redundant rows of a given parity check matrix was discussed 
in [13] [4]. In their methods, the redundant rows are effi- 
ciently found based on a short cycle of a Tanner graph of 
the parity check matrix. Their results indicate that addition 
of redundant rows is a promising technique to improve LP 
decoding performance and they suggest that further studies on 
this subject are meaningful to pursue more systematic ways to 
find appropriate redundant rows that achieves better trade-offs 
between decoding performance and complexity. 



II. Preliminaries 

In this section, notations and definitions required throughout 
the paper are introduced. 

A. Notations and definitions 

Let H he a binary m x n matrix where n > m > I. The 
binary linear code defined by H, C{H), is defined by 

C{H) = {xe F^ : xH^ = 0™} , (1) 

where F2 is the Galois field with two elements {0, 1} and 0™ 
is the zero vector of length to. In the present paper, the bold 
face letter, like x, denotes a row vector. The elements of a 
vector is expressed by corresponding normal face letter with 
subscript; e.g., x — {xi, X2, • ■ • , Xn)- The following definition 
of the fundamental polytope is due to Feldman [3]. 

Definition 1: (Fundamental polytope) Assume that t is a 
binary row vector of length n. Let 



X{t) = {SC Supp(t) : |5| is odd}, 



(2) 



where Supp(t) denotes the support set of the vector t defined 
by 

Supp{t) = {ie{l,...,n}:U^O}. (3) 

The single parity polytope of t is the polytope defined by 

W(t) = {/e[0,l]":V5eX(t), 

E /j- + E (1 - fi) ^ |Supp(*)l - 1}, (4) 



jes 



je(Supp(t)\s) 



where [a,b] = {v G M. : a < v < b}. The symbol M denotes 
the set of real numbers. For a given binary mxn parity check 
matrix H, the fundamental polytope V{H) is defined by 



V{H) = f]U{h,), 



(5) 



where hi — {hn, hi2, . . . , hin) is the i-th row vector of H. q 
The convex hull of C{H) is the intersection of all convex 
sets including C(i7r] It is known that the fundamental 
polytope V{H) contains the convex hull of C{H) as a subset 
and that V{H) n {0, 1}" = C{H) holds. 

Let M denote the number of the linear constraints (inequal- 
ities) that forms V{H). Assume that these linear constraints 
are numbered from 1 to M and that the fc-th linear constraint 
has the form; 



Oiklfl H 1- Oiknfn < /3fc, 



(6) 



for fc e {1, . . . , A/}. We call the fc-th constraint Const^. The 
hyper plane corresponding to Constfc is given by 

Tk ^ {/ G M" : afci/i + ■ • • + aunln - h}. (7) 

and the half space satisfying Constfc is defined by 

nk = {fe M" : afei/i + • • • + aknfn < Pk}- (8) 

'in this context, we implicitly assume that C{H) is embedded in K". 



The fundamental polytope 1^(11) is thus the intersection of 
the half spaces such that 



r{H) = Til n W2 n • • • n 7^ 



M- 



(9) 



Let Fact(-ff) be the set of indices of the active constraints 
defined by 



F,,tiH) = {ke{l,...,M}:0" GTk}- 



(10) 



The active constraints are the linear constraints whose hyper 
plane contains the origin 0". In a similar manner, we define 
Finact(-H^), which is the set of indices of the inactive con- 
straints, by 



tiH) = {ke{l,...,M}:0" (^Tk}. 



(11) 



The fundamental cone JC{H) is the cone defined by the active 
constraints: 

K.{H) = {fe [0, 1]" : Vfc e Fact(i/), / e Hk}- (12) 
III. Cutting plane method 

In this section, the main idea of the cutting plane method 
based on redundant rows will be introduced and then an 
application to Hamming code is shown as an example. 

A. Fractional distance 

The fractional distance dfindH) is the £1 -distance between 
a codeword and the nearest vertex of 1^(11) [3]. 

Definition 2: (Fractional distance) Let V{X) be the set of 
all vertices of a polytope X. For a given binary mxn parity 
check matrix H, the fractional distance of H is defined by 



-^frac 



(H) 



min > \xi 

x<£C{H) 

feviv(H)) ' 



i=l 



M- 



(13) 



D 



The importance of the fractional distance is stated in the 
following lemma due to Feldman [3]. 

Lemma 1: Assume that the channel is BSC. Let e be the 
number of the bit flip errors. If 



dfiixc{H) 



1 



>e, 



(14) 



holds, then all the bit flip errors can be corrected by the LP 

decoder 

(Proof) The proof is given in [3]. q 

Based on the geometrical uniformity of the fundamental 

polytope (called C-symmetry in [3]), it has been proved that 

dfrac{H) is the ^i-weight of the minimum weight vertex of 

V{H) except for the origin, which is expressed by 



dft.ac(i/)= min > /,. (15) 

/ev-P H \0" ^ 

?.— 1 

Let T{H) be a set of the minimum weight vertices of V{H)\ 
T{H) = lpe V{V{H)) : Ep* = rffrac(i?) \ , (16) 



and let (ijjji„ be the minimum distance of C{H). Since 
any codeword of C{H) is a vertex of V{H), it is obvious 
that the inequality d{inc{H) < d^ain holds for any H. The 
fractional distance (ifiac(^) depends on the representation of 
a given binary linear code (i.e., the parity check matrix) and 
there are a number of parity check matrices that define the 
same binary linear code. This means that the parity check 
matrices of a binary linear code can be ranked in terms of its 
fractional distance. It is hoped to find a better parity check 
which achieves larger factional distance for a given binary 
linear code because such a parity check matrix may improve 
the performance/complexity trade-off of LP decoding for the 
target code. 

B. Cutting polytope 

In this subsection, we will define the cutting polytope based 
on the redundant rows. The next definition gives the definition 
of the redundant row. 

Definition 3: (Redundant row) Let H he a binary rn x n 
parity check matrix of the target code C. A redundant row h 
is a linear combination of the row vectors of H such that 



h = aihi + 02/12 



7, *''n 



(17) 



where Ci e F2(i e {!,. ..,m}). q 

The single parity polytope of a redundant row h includes all 
the codewords of C{H) because any codeword x e C{H) 
satisfies xh^ = [12]. 

The next lemma asserts a cutting property of a single parity 
polytope satisfying a certain condition. 

Lemma 2: Let i/ be a binary m x n parity check matrix. 
Assume that a point p = (pi, . . . ,p„) G M" and i G Fg are 
given. If 



3j e Supp(<), 



iesupp(t)\{i} 



(18) 



holds, then p ^ U{t) holds. 

(Proof) Let j* be the index satisfying 

Pr > H Pi- 

ieSupp{t)\{j'} 



(19) 



The above inequality is equivalent to the following inequality: 

(l-pO>|Supp(<)|-l. (20) 



Y. Pi 
je{j*} 



E 

iGSupp(t)\{i-} 



From the definition of U{t), it is evident that p ^ U{t). \^ 
The cutting polytope defined below is used to cut a 

fractional vertex in the proposed method for improving the 

fractional distance. 

Definition 4: (Cutting polytope) Assume that p ~ 

(pi,P2, ■ ■ ■ iPn) € T^{H). If a redundant row h with the form: 



satisfies 



h = aihi + 02/12 
Pj > 



i.'^n 



E Pi, 

ie(Supp(h)\{j}) 



(21) 
(22) 



where j — argmaXjggupp(^-)Pj, then U{h) is called a cutting 
polytope of p. Q 

The next theorem introduce a tighter relaxation of the con- 
vex hull of C{H) which may improve the fractional distance. 

Theorem 1: Let p — {pi,p2, . . . ,Pn) G 'PiH) and U{h) 
be a cutting polytope of p. The following relations hold: 



and 



ciH)cr{H)nuih) 



piV{H)r\U{h). 



(23) 



(24) 



(Proof) The claim of the theorem is directly derived from 
C{H) C U{t) and Lemma [I] g 

The theorem implies that the cutting polytope of p also 
contains C{H) but excludes p. In other words, the intersection 
V{H) n U{h) is a tighter relaxation of the convex hull of 
C{H) compared with V{H). Let H' be the parity check matrix 
obtained by stacking H and h, namely. 



H' 



H 
h 



(25) 



Note that 'P(-ff') = V{H) D U{h) has also geometrical 
uniformity because V{H') is a fundamental polytope. This 
means that the cutting polytope cuts not only p but also 
some non-codeword vertices ofV{H) which are geometrically 
equivalent to p. The fractional distance of H' , dii-a,c{H'), thus 
can be larger than the fractional distance (ifrac(-H^) because 
the point p E r(iJ) is excluded from the new fundamental 
polytope V{H'). Furthermore, the multiplicity (i.e., size of 
V{H)) can be reduced as well by eliminating a fractional 
vertex with the minimum £i-weight. Figure [T] illustrates the 
idea of the cutting polytope. 



Codeword 




^>^ Cutting 
polytope 
U{h) 



Fig. 1 . The idea of tlie cutting polytope 



C. Cutting plane method: an example 

In this subsection, we will examine the idea of the cutting 
plane method described in the previous subsection thorough a 
concrete example. Let i/ be a parity check matrix of (7, 4, 3) 



Hamming code: 






/ 1 1 1 1 


H = 


110 10 10 




I 1 1 1 1 



(26) 



In this case, we have the index sets: 

Supp(/ii) = {l,3,4,5}, 
Supp(/i2) = {l,2,4,6}, 
Supp(/i3) = {2,3,4,7} 

and 

X(/ii)^{{l},{3},{4},{5}, 

{3,4,5},{1,4,5},{1,3,5},{1,3,4}}, 
X(/^2) = {{1},{2},{4},{6}, 

{2,4,6},{1,4,6},{1,2,6},{1,2,4}}, 
X(/i3) = {{2},{3},{4},{7}, 

{3, 4, 7}, {2, 4, 7}, {2, 3, 7}, {2, 3, 4}}. 

The fundamental poly tope of H is the set of points in [0, 1]^ 
satisfying 






E 

je(Supp(hi)\s) 



(1 - lo) < 3 



(27) 



for any i E {1,2,3} and any S E X{hi). From some 
computations (details of computation are described in the next 
section), we can obtain the set of £i -minimum weight vertices 
ofV{H)({.t.,T{H)): 



0,^,^,0,0,0 



^,o,y,o,o,o 



y,o,|o,o,o 



Therefore, in this case, di^s,c{H) is equal to 2. 

Assume that we choose h = (1, 0, 1, 0, 0, 1, 1) as a redun- 
dant row that is the sum of the second and third rows of H. 
Let p = (0, 2/3, 2/3, 2/3, 0, 0, 0) E T{H). It is easy to check 
that 

"^ (28) 



P3 



> 



;esupp(h)\{3} 



Pi 



holds where Supp(/i) = {1,3,6,7}. This means that U{h) 
is a cutting polytope of p. By stacking H and h, we get a 
new parity check matrix H' whose fundamental polytope does 
not contain p as its vertex. In a similar manner, continuing the 
above process (appending redundant rows to H for cutting the 
vertices in T{H)), we eventually obtain a parity check matrix 

H*: 

/ 1 1 1 1 \ 

110 10 10 
1110 1 



H* = 



1 
1 


Vo 



1 
1 
1 

0/ 



(29) 



The fractional distance of H* is equal to 3 which is strictly 
larger than the fractional distance of H {dhadH) — 2). The 
details of a way to find an appropriate redundant rows will 



be discussed in the subsequent sections. It has been observed 
that the vectors in T{H*) are integral; namely, they are the 
minimum weight codewords of Hamming code. 

D. A greedy algorithm for cutting plane method 

The previous example on Hamming code suggests that 
iterative use of the cutting plane method for a given parity 
check matrix may yield a parity check matrix with redundant 
rows which is better than the original one in terms of the 
fractional distance. The following greedy algorithm, called 
greedy cutting plane algorithm, is naturally obtained from the 
above observation. 



Greedy cutting plane algorithm 




Step 


1 Evaluate T{H). 




Step 


2Pickupper(ff). 






Step 


3 Find a redundant row 


h of H which 


gives a 




cutting polytope of p. 


If such h does not exist. 




exit the procedure. 






Step 


4 Update H by 


( H\ 






H := 


\ h )■ 


(30) 


Step 


5 Return to Step 1 . 







Details on the process of finding a redundant row h is 
shown in the next section. The most time consuming parts 
of the above algorithm are evaluation of T{H) and search for 
a redundant row. In the next section, we will discuss efficient 
implementations for these parts that are indispensable to deal 
with codes of long length. 

IV. Efhcient implementation 
A. Efficient computation of dhii.c{H) and T(H) 

As described in the previous section, evaluation of dfrac(^) 
and T{H) is required for finding a redundant parity check ma- 
trix with better fractional distance. An algorithm for computing 
diTac{H) has been proposed by Feldman [3]. 

Firstly, we review the Feldman's method. For any k E 
{1,...,M}, Let dk be 



dk = minimize ^ fi s.t. / E {V{H) D Tk) 



(31) 



where this LP problem is denoted by LPk- Thus, dk can be 
considered as the ^i -weight of the minimum weight vertex on 
the facet Tk of V{H). Since there exists at least one facet of 
V{H) which includes p for any vector p in T{H), it is evident 
that 



.(H) 



M 

mmSk 

fe=i 



(32) 



holds where Sk = dk if dk > 0; otherwise Sk = oo. The 
LP problems LPk can be efficiently solved with an LP solver 
based on the simplex algorithm or the interior point algorithm. 
From the solution of these LP problems, we can obtain the 



fractional distance of i/ by ( 32 1 



The number of constraints related to a fundamental polytope 
defined by (J5]l is an exponential function of the row weight of 
H. Thus the number of executions of the LP solver rapidly in- 
creases as the row weight of H grows. Another formulation of 
the fundamental polytope proposed by Yang et al.[15] can be 
used to reduce the number of constraints. In their formulation, 
high weight rows of H are divided to some low weight rows 
by introducing auxiliary variables. Although their method is 
effective for evaluating the fractional distance as well, we here 
propose another efficient method for evaluating dfrac(^) in 
this section. In our method, a fundamental polytope is relaxed 
to a fundamental cone. This method can be combined with 
Yang et al.'s formulation. 

We need to prepare a relaxed version of the fractional 
distance before discussing another expression of the fractional 
distance. For k e Fi,iact(-ff). let f^^ be 



^Max) ^ ^^— ^^-^c Y^ f^ s_t^ f (= (^^H) n Tk). (33) 



This relaxed LP problem is denoted by LP^^'^ ^'^ . The relaxed 



fractional distance d^ 

.(relax) 
"frac 



(relax) 



is defined by 



(H) 



keFi„^ct(H) 



,(rclax) 
^k 



(34) 



The next theorem states a useful equivalence relation. 
Theorem 2: For a given H, the following equality holds: 



dt^r\H)^d,,uH). 



(35) 

(Proof) See appendix. D 

A merit of Theorem E is that the evaluation of rffraJ*'^ (^) 
takes less computational complexity than that of the evaluation 
of offrac(^) using the Feldman's method. The reduction on 
the computational complexity comes from the following two 
reasons. 



(relax) 



One reason is that the feasible region of LP^ ' is 
based on the fundamental cone K.{H) (instead of V{H)) 
which can be expressed with fewer linear constraints than the 
fundamental polytope. In the case of a regular LDPC code 
with row weight w^, the number of linear constraints required 
to define ipj-'^'^) 



is mwr + L On the other hand, to express 
the fundamental polytope required in £Pfe, ni2'^^^^ linear 
constraints are needed. 

Another reason of the reduction on complexity is that 
fewer executions of the LP solver are required for evaluating 



LP, 



(relax) 



because we can focus only on the inactive linear 
constraints in the case of LPj^^ . 

B. Search for redundant rows 

A straightforward way to obtain a redundant row h that 
gives a cutting polytope of a given point p is the exhaustive 
search. Namely, each redundant row is checked whether it sat- 



isfies the condition ( 18 1 or not. However, this naive approach is 
prohibitively slow even for codes of moderate length because 
there are 2™ redundant rows. We thus need a remedy to narrow 



the search space. The following theorem gives the basis of the 
reduction on the search space. 

Theorem 3: Assume that there exists a cutting polytope of 
p E r(iJ). The polytope U{h*) is such a cutting polytope of p 
and the redundant row h* is given by h* — X]"=i o*ihi. Then, 
U{h) is also a cutting polytope of p where h = X]"=i ^t^i 
and 

a* ieQ 
i4Q. 



(36) 



The index set Q is defined by 

g = {z e {1, . . . , m} : 3j e Supp(p), /i,j ^ 0} . (37) 
(Proof) Let j* e Supp(ft,*) n Supp(p) be the index satisfying 

/esupp(h*)\{j*} 

Since pj = if j ^ Supp(p), the above condition is equivalent 
to 

ie(Supp(h-)nSupp(p))\{j-} 
From the definition of h and Q, we have 

Supp(/i*) n Supp(p) = Supp(/i) n Supp(p). (40) 

This equality leads to the inequality 

Pj' > Y Pi 

ie(Supp(h)nSupp(p))\{i-} 

Y. pi- (41) 

iesupp(h)\{j-} 

The above inequality implies that U{h) is a cutting polytope 
of p. □ 

The significance of the above theorem is that we can fix 
fli = for i ^ Q in a search process without loss of the 
chance to find a redundant row generating a cutting polytope. 
Therefore, computational complexity to find a redundant row 
can be reduced by using this property. Let V be a sub-matrix 
of H composed from the columns of H corresponding to the 
support of p. The index set Q consists of the row indices of 
non-zero rows of V. Thus, in the case of LDPC codes, the 
size of Q is expected to be small when the size of Supp(p) 
is small because of sparseness of the parity check matrices. 
In such a case, the search space of the redundant rows are 
limited in a reasonable size. For example, in the case of the 
LDPC code "96.33.964" [9](n = 96,to = 48, row weight 6, 
column weight 3), the size of Supp(p) is 7{p E r(iJ)) and 
the size of Q is 8. 

Let H'^ be the |(5| x n sub-matrix of H composed from a 
row vectors whose indices are included in Q. From Theorem 
3, we can limit the search space to the linear combinations 
of rows of H'^. In the following, we will present an efficient 
search algorithm for a redundant row. Let x = (ti, . . . ,t„) 
denote an indices vector that satisfies p^^ > p^-o > ■ ■ ■ > Pt„ ■ 



Redundant row search algorithm 

Step 1 Construct H'^. 

Step 2 Permute columns of H'^ to the following form: 

H'^U = ( Vr^ Vrs ■■■ Vr,^ V^, ) , 

where 11 denotes a column permutation matrix 

and Vj,j G {1, . . . , n} denotes the j-th column 

vector of H'^ . 
Step 3 Convert H'^H into U of row echelon form by 

applying elementary row operations. 
Step 4 Let 

i* = argminji G {1, 

where Ui = {un, . 
vector of UI1~^. 
Step 5 Output Ui* . 



,\Q\} : u, satisfies Q}, 
, Uin) denotes the i-th row 



The idea of the redundant row search algorithm is based on 
the fact that Ui (a candidate of desirable redundant rows) tends 



to satisfy condition (22i. This can be explained as follows. 



Assume that Ui^n = 1- From the definition of row echelon 
form, 

holds where K is an integer larger than or equal to i — 1. This 
means that 

|Supp(M,) n Supp(p)| < |Supp(p)| - i + 1 (42) 

holds for i e {1, . . . , \Q\}. Let r^., be 



E 



pi- 



(43) 



ie(Supp(ui)nSupp(p))\{ri} 



From the inequality (42 1, it is evident that {?7i,7?27 ■ • •} is 
a decreasing sequence. We thus can expect that condition 
(22 1, i.e., pr^ > rji, eventually holds as i grows^ It may 
be reasonable to choose the smallest index i satisfying (22i 



because such Ui would be sparser in the case of a low density 
matrix. A sparse redundant row is advantageous since it is able 
to cut other fractional vertices with small weight. 

V. Results 

A. Application to Hamming, Golay and LDPC codes 

In this subsection, we applied the cutting plane method 
to Hamming code (n — 7,rn — 4), Golay code (n = 
24, m = 12), and a short regular LDPC code "204.33.484" 
[9] (n = 204, m — 102, row weight 6, column weight 3). 
We here use a parity check matrix of Golay code described in 
[12]. Let dfrac be the fractional distance of the original parity 
check matrices and df/*" be the fractional distance of parity 
check matrices generated by the cutting plane method. Let Nd 
be the number of the appended redundant rows. The results 
are shown in Table lit) For example, addition of 100 rows to 
the original parity check matrix of Golay code increases the 

-Note that we cannot state that Ui satisfying 122} always exists because 
this argument is based on the assumption Ui ^-^ =1. 



TABLE I 

Fractional distances of redundant parity check matrices 
obtained by the cutting plane method 



Code 


n 


m 


Nd 


dfrac 


patter 


f^niin 


Hamming 


1 


3 


4 


2.000 


3.000 


3 


Golay 


24 


12 


40 


2.625 


3.429 


8 


Golay 


24 


12 


100 


2.625 


3.895 


8 


LDPC 


204 


102 


9 


5.526 


5.646 


unknown 



o 



k^\ 


original -<»- 

+40rows -<•- 

+100rows — •— 




\ \ \. 


\ ^- 


\ X 


^v^ 


X^ 


^^ 




■* 



Crossover Probability on BSC 

Fig. 2. Comparison on block error probabilities of parity check matrices of 
Golay code (original 24 X 12 matrix, redundant 24 X 52 matrix, redundant 
24 X 112 matrix ) 



fractional distance from 2.625 to 3.895. It is expected that 
these matrices constructed by the cutting plane method shows 
better LP decoding performance compared with the original 
matrices because 



'^ftac'^ is greater than dfrac for all the cases. 



B. Decoding performances 

In this subsection, we will present decoding performance of 
redundant parity check matrices obtained by the cutting plane 
method. We here assume BSC as a target channel and LP 
decoding [3] as a decoding algorithm used in a receiver Figure 
[2] shows block error probabilities of the 24 x 12 parity check 
matrix of Golay code given in [12] (labeled "original"), the 
24 X 52 parity check matrix with 40 redundant rows (labeled 
"+40rows"), and the 24 x 112 parity check matrix with 100 
redundant rows (labeled "+100rows"). 

From Figure l2j we can see that the block error probability of 
the redundant matrix (+100rows) is approximately two order 
of magnitude lower than that of the original matrix when a 
crossover probability is 10^^. Figureplshows block error prob- 
abilities of the original parity check matrix of "204.33.484" 
[9] and that the parity check matrix with 9 redundant rows 
(labeled "+9rows"). From Figure l3] it can be observed that 
the slope of the error curve of "+9rows" is steeper than that 
of "original". 

VI. Conclusion 

In this paper, the cutting plane method based on redundant 
rows of a parity check matrix for improving the fractional 
distance has been presented. In order to reduce the search 
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Fig. 3. Comparison of bloclc error probability of parity check matrices of 
the LDPC code (n = 204, m = 102) 



complexity to find an appropriate redundant row, we intro- 
duced an efficient technique to compute the fractional distance 
and proved that the limited search space indicated in Theorem 
|3] is sufficient to find a redundant row generating a cutting 
polytope. Some numerical results obtained so far are encour- 
aging. The redundant parity check matrices constructed by the 
cutting plane method have larger fractional distance than that 
of the original matrices. The simulation results support that 
improvement on the fractional distance actually leads to better 
decoding performance under LP decoding. 

Appendix 

Proof 1 (Theorem 1): For proving Theorem 1, we will 
prove the following two inequalities: 

dtfr\H) > d,uH). 



(44) 
(45) 



7(rclax) 



We first assume df^^^^ (H) > dfrac(-ff) for proving in- 
equality ( |44l ) by contradiction. Assume that the index i E 
{!,..., il/} is given by i e argmin^.^]^ (5^. Let p be the 
solution of LPi. This means that p € r(i?) and p g {V{H) n 
Ti) where Vi^H) n Ti is the feasible set of LPi. Note that 
p is not the origin 0" because 6k = oo holds when p = 0". 
Thus, Const j is an inactive constraint, namely i E Finnct{H). 
Since the fundamental corn K.{H) contains the fundamental 
polytope 7^(i?) as a subset, the relation p E {JC{H) n JFj). 
holds as well. The intersection JC{H) n Ti is the feasible 
set of LP^'^^'. Let the solution of Lpj-^'^ ^^' be p* . Since 

(relax) 

and the feasible set of LP^ includes that of LPi, the 

£i -weight of p* is smaller than or equal to the ^i -weight of 
p. This implies d^^'^^ {H) < dfi.ac(-ff) but it contradicts the 
assumption. The proof of the inequality ( |44| is completed. 

(relax) 



the objective function of LPi is identical to that of LP^ 



We next assume d\^^^ {H) < 



dha,c{H) for prov- 
ing inequality ( 45 1 by contradiction. Assume that i E 

argminfcgFi„„t(if)"'fe 



^(relax) j^^^ ^(relax) ^,g jjjg golution of 



p(reiax) ^ {K.{H) H T,) holds. In the following, we will 
discuss the two cases: (i) p('<=ia'') e V{H), (ii) p('<=ia'=) ^ 
r{H). We start from case (i). If p('<=iax) g p(^) ^olds, 
then p(''<=iax) g (7'(7J) n Ti) holds because of the relation 
P{H) C IC{H). Due to almost the same argument, we obtain 

^(relax) 

then move to case (ii). If p(''<=iaxj 
Const; satisfying 



dfl'^'^^'{H) = d{,-iic{H) which contradicts the assumption. We 



^ 'P{H) holds, there exists 



(relax) 



p— iUu 0" e H/, IE Fi„aet(i?). (46) 

This is because pi^^^^^) e JC{H) but ^(^1^^^) ^ V{H). Let 
Seg(p*^'''''^''') be the line segment between p(''eiax) ^^^ ^^ 
origin: 

Seg(p("'"")) = {pei?":0<t<l,p = ip('°'"'^)} . (47) 

The line segment Seg(p'^'''''^''^) passes through Ti. Thus, there 
exists the point pi E {Ti n Seg(p(™''"'))). Note that the line 
segment Seg(p'^''°'^'^)) is totally included in IC{H) because 
p(reiax) ^ ^^y rj,j^jg jg^j^ j^ ^^ ^ {IC{H) C] Ti) and it 

means that pi is included in the feasible set of LP; 



(relax) 



The 



£i -weight of pi is smaller than that of p(.'<=iaxj ^\^^ implies 
that the ^i-weight of p; is smaller than d^°^^ {H). However, 

(relax) 
^frac 



it contradicts the definition of de^'^T'^'{H). The proof of the 

n 



^p(ieax)^ From the definition of l^p^^'^^' jt is evident that 



inequality (45i is completed. 
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